Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1–13, 15, 17–20, 24, 26–28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.


Trait records
Weekly body weight data of 805 quail (399 females and 406 males) from the parental, F 1 and F 2 generations were collected from hatching until 16 weeks of age (WK 0-16).A total number of 3990 eggs obtained from 399 females were used to assess egg-related traits.Egg-related traits were measured at two different egg production stages: at the beginning of the egg production stage (first stage) and at 12 weeks of age (second stage).The two stages are indicated with subscript letters ' 1 ' and ' 2 ' in the abbreviations of traits.The first five eggs from each female were evaluated at each developmental stage.External and internal egg-related traits, including egg weight (EW), egg long axis (ELA), egg short axis (ESA), eggshell strength (ESS), eggshell weight (ESW), egg equator thickness (EET), yolk weight (YW), yolk diameter (YD), yolk colour, lightness (L* value) (YC-L*), redness (a*

Association analysis
Association analysis was performed using general linear model (GLM) and mixed linear model (MLM) approaches 63 implemented by TASSEL v. 5 64 .GLM evaluation incorporated population structure (Principal component analysis, (PCA)) and MLM used both population structure (PCA) and kinship matrix in the association analysis.The PCA and covariance analysis were performed using the method reported by Ref. 65 to describe the population structure.A kinship matrix was produced using the Kinship tool with the Scald_IBS method to avoid spurious associations due to relatedness and population structure.In the "Scald_IBS" method, genotypes were encoded based on the count of the alleles at that locus.By utilizing this encoded data, the kinship matrix was estimated.The efficient mixed-model association (EMMA) algorithm 66 and population parameter previously determined (P3D) variance component estimation were implemented to decrease the computing time for a large dataset 67 .Statistically significant associations with p-values < P threshold = 0.01/N, where N is the total number of SNP markers, were identified, and a standard correction was performed by applying a Bonferroni procedure at a minimum value of 3.88806E −36 for GLM and 7.03385E −22 for MLM 68 .The coefficient of determination (R 2 ) was subsequently calculated to quantify the extent to which genetic markers explained the phenotypic variance.The significant value and the marker effect for each SNP were exported, and Manhattan as well as quantile-quantile (QQ) plots were generated in the R Project for Statistical Analysis v. 3.6.1 69 using the qqman package 70 .The data were subjected to a correction procedure to account for the influence of fixed effects, specifically birth date, sex, and dame.This adjustment was performed using the least squares method in R-Project for Statistical Analysis v. 3.6.1 69 .

Ethical approval
Animal care, experimental protocols, and blood sample collections were approved and conducted in accordance with the Rules on Experimental Animals and Animal Experiments at Hiroshima University, Graduate School of Integrated Sciences for Life, Laboratory of Animal Breeding and Genetics (Approval No. C20-15) and the protocol described in the Guidelines for Proper Conduct of Animal Experiments, Science Council, Japan https:// www.scj.go.jp/ ja/ info/ kohyo/ pdf/ kohyo-20-k16-2e.pdf.Also, we confirm that all methods have been conducted in adherence to the ARRIVE guidelines (https:// arriv eguid elines.org).

Results
Filial generations obtained from reciprocal matings are presented separately.The body weights of female and male birds in the LS and NS strains and their F 1 and F 2 generations are shown in Tables 1 and 2, respectively.In females, LS F 1 (LS♀ ✕ NS♂) quail showed the highest body weight at hatching (8.19 and 8.54 g, respectively).The lowest body weight was detected in NS quail from hatching to 16 weeks of age.No difference was found between F 2 generations in WK 0-2.From 4 to 16 weeks of age, the LS strain presented the highest body weight compared to that of the other groups.In males, the LS F 1 (LS♀ ✕ NS♂) quail showed the highest body weight at birth.However, no difference was observed between the weights of the other quail groups, except for NS.Moreover, no significant differences were observed between the F 1 generations in WK 1-16.The NS quail showed the lowest body weight during the entire period analysed.The LS strain exhibited the highest body weight from 4 to 16 weeks of age compared to those observed at all generations.
Table 3 shows the means and standard errors for egg-related traits in the LS and NS strains of Japanese quail and their F 1 and F 2 hybrids in the first and second egg-laying stages.The LS strain EW, ESA, ESW, and AW exhibited the highest value among all groups in both egg-laying stages.NS females showed the highest values for the YC-a* 1 and YC-b* 2 traits across generations.F 2 (LS♂ ✕ NS♀) females showed the lowest value in YC-a* 1 (2.64) compared to the those observed in the other groups.No difference was observed in EET 1 between the parental and filial generations.No significant differences were observed in ELA 1 between NS and both filial generation females.Both F 2 groups presented the lowest ESS 1,2 among all groups, whereas no differences were observed between the parental and F 1 quail in the first and second stage of egg production.F 2 (LS♂ ✕ NS♀) birds started  Illumina HiSeq 2500 produced 733,233,057 RAD-seq reads for the treated samples.A total of 25,631 SNPs that passed the quality control were used for the GLM and MLM analyses.This subset of SNP genotypes was used to produce a kinship matrix for MLM analysis.The same subset of SNP markers was also used to perform the PCA in this study.General and mixed linear models produced 1,941 and 2,986 associations, respectively, after applying statistically significant loci with Bonferroni corrections.Statistically significant associations with p-values lower than the threshold were excluded, and 383 SNPs for 1,151 associations (Supplementary File S1) and 734 SNPs for 1,442 associations (Supplementary File S2) were identified to be significantly associated with production traits in the GLM and MLM approaches, respectively.GLM identified SNPs that were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying all analysed traits except EW 1 , AW 1 , and YC-b* 2 (Table 4).The maximum number of SNPs and associations was identified on the Z chromosome (122 and 611, respectively), followed by chromosome 1 (61 and 95, respectively) and chromosome 4 (52 and 92, respectively), while the minimum number of SNPs and associations (1 and 1, respectively) was found on chromosomes 13, 20, and 28 underlying WK 0, ESS 1 , and WK 3, respectively.A total of 1065 and 86 associations were detected for body weight and egg-related traits, respectively (Table 4).No associations for body weight traits were identified on chromosomes 12 (5,472,506 bp) and 20 (11,476,337 bp).Additionally, no associations for egg-related traits were detected on chromosomes 13, 24, and 26-28.The percentage of phenotypic variance explained by marker (R 2 ) detected by the associations ranged between 0.003 and 0.097 using the GLM approach.Figures 1 and 2 represent the Manhattan and QQ plots for ESS 2 , respectively, wherein 14 associations were found on chromosomes 1, 2, 4, 5, 11, and Z.The Manhattan and QQ plots for the remaining traits from the GLM are shown in Supplementary File S4.
As shown in Table 5, MLM recognised SNPs that were positioned on all chromosomes except for chromosome 16 and associated with all traits except for ELA 2 .The maximum number of associations was detected on the Z chromosome (N = 516), followed by chromosomes 1 (N = 208), 3 (N = 105), and 2 (N = 96), while the minimum number of associations was observed on chromosome 25 (N = 1) for YC-L* 2 and chromosome 28 (N = 1) for egg production rate.However, a maximum number of 149, 97, 76, and 69 SNPs were distributed on chromosomes Table 3. Means and standard errors for egg-related traits in LS and NS strains of Japanese quail and their F 1 and F 2 hybrids.a-e Means with different superscript letters show the significantly different in all individuals (Tukey's HSD test, P < 0.05).1,2 Subscript letters are first and second egg laying stages.The data of LS, NS, F 1 (LS♂ ✕ NS♀), and F 2 (LS♂ ✕ NS♀) are the same as those in our previous paper 17 .and 4, respectively.For WK 1 trait, four associations were detected on chromosomes 1 (N = 2), 3 (N = 1), and Z (N = 1).MLM plots for each Manhattan and QQ (Supplementary File S5) were plotted for the associated body weight and egg-related traits.

Discussion
The evaluation of 44 phenotypic traits from 567 birds in the F 2 revealed 383 SNPs with 1151 associations and 734 SNPs with 1442 associations for GLM and MLM approaches, respectively.Our identified SNPs were associated with all targeted traits except EW 1 , AW 1 , and YC-b* 2 in the GLM and ELA 2 in the MLM analysis.To the best of our knowledge, no such association analysis of production traits in Japanese quail have been reported using GLM and MLM approaches with RAD-seq data.However, several studies on QTL analysis of production traits in Japanese quail have been reported.Reference 20 recognised QTLs underlying body weight at hatching and at four weeks of age on chromosome 1 at 12-13 cM.Here, we identified associations on the same chromosome at 11.9 cM for GLM and 12.38 cM for MLM corresponding for WK 0 and TLE traits, respectively.Reference 21 reported QTL associated with body weight from hatching to three weeks of age on chromosome 3 between 27 and 30 cM.The results of the present study revealed trait associations for SNPs located between 28.73 and 37.74 cM on chromosome 3 in both GLM and MLM analyses.Moreover, using GLM and MLM approaches, we showed that associations underlying body weight traits were mostly located on chromosomes 1-10 and Z. Associations affecting hatching weight were identified on chromosomes 1-9, 13, 17, 27, and Z for GLM, and 1, 3, 7, 9, 13, and 27 for MLM approaches.Detected SNPs located at 20.16 cM on chromosome 7 were associated with hatching weight in both GLM and MLM analyses, and these results were consistent with the previously identified location of QTLs associated with this trait on chromosome 5 (19-20 cM) 22 .The detected QTL affecting liveweight measurement at week 5 has been shown to be located on chromosome 1 in a study investigating the genetic mapping of QTLs affecting body weight 18 , which is in agreement with the SNPs associated with body weight and located on chromosome 1 in our GLM and MLM analyses.Reference 3 identified a QTL for body weight at 4 weeks of age at the centromere of chromosome 2 which is consistent with the locus detected for this trait on chromosome 2 in both association analyses performed in the present study.The study detected three QTLs associated with body weight at 1, 4, and 6 weeks of age, located in the initial region of chromosome 2 between 0 and 15 cM.Reference 2 reported putative pleiotropic loci on chromosome 3 (52.6-56.7 Mb) affecting both weight and egg traits, where a QTL for egg weight co-localised with a QTL for body weight at 65 days of age.Interestingly, our results indicated an association in the same chromosome (56.08Mb) controlling WK 1 in both the GLM and MLM analyses.This candidate region could underlie the genetic correlation already observed in quail between body weight and egg traits [71][72][73] .QTLs for early and late growth stages (5 and 70 weeks of age) have been previously detected on chromosome 1, positioned 18-19 cM 24 .In our MLM analysis, YW associations were identified on chromosome 1 between 17.68 and 21.16 cM.It is likely that the QTLs for BW and YW located on chromosome 1 represent a single gene.In addition to QTLs for body weight traits, Ref. 24 identified QTLs for eggshell weight on chromosomes 1 (191 cM), 5 (12 cM), and 20 (21 cM); egg weight, egg number until the age of 69 weeks, and age of first egg on chromosome 6 positioned at 0, 32, and 34 cM, respectively.The results of the present study revealed associations for YC-L* 2 and YC-a* 2 on chromosome 1 (190.34cM) as well as TLE and EPR on chromosome 6 (31.8 cM) in GLM.Alternatively, we identified EET 2 QTLs on chromosome 5 (12.02 cM), YW 2 on chromosome 6 (0 cM), and YW 2 as well as YD 2 on chromosome 20 (20.35 cM) in the MLM.These findings are consistent with the results presented by Ref. 24 .In turn, Ref. 3 recognised QTLs for the number of eggs laid and the egg production rate on chromosome 1 between 36 and 42 cM.The differences among the results of these studies may be explained by the different durations of egg production.In our previous report, we identified a QTL for growth-related traits on chromosome 1 using the same Japanese quail strains 16 .Similarly, in the present study, two associations were found on the same locus (150,370,977 bp) on chromosome 1 for BW 4 and 5.In addition, we detected QTLs for ESA 1 and YD 1 on the Z chromosome between 22,757,726 and 31,279,475 bp 17 .Associations underlying WK 9 in GLM, EET 1 in MLM, and WK 7-16 in both approaches were identified within the same chromosome and at the same positions.This was expected as positive correlations have been observed between the related egg trait and body weight 73 .Only chromosomal regions were explored in the above-mentioned studies on the Japanese quails; therefore, future studies focusing on the association between SNPs and candidate genes would be of great value for improving QTL resolution in Galliformes.GWASs are widely used to distinguish SNPs associated with production traits in chickens 74,75 .In the first study of its kind on Japanese quail, we used association analysis with GLM and MLM, based on RAD-seq SNPs, to examine production traits like body weight and egg-related traits.Comparative genomic studies of Japanese quail and chickens based on cytogenetics 53 , orthologous genes 76 , and linkage analysis 52 have shown a high rate of synteny-conserved karyotypes and genomic structure between these species.Therefore, the GWAS results from chickens were used to identify potential candidate genes in the present study.However, due to differences in growth conditions between chickens and quail, body weight at the early growth stage, sexual maturity age, and maximum growth rate were considered in both species.
The detection of associations related to body weight during juvenile age is likely due to the strong correlation between body weight traits during early growth stages.It is anticipated that numerous genes control body weight during this period because growth is a highly intricate trait influenced by multiple loci that affect appetite, feed intake, body composition, nutrient utilization, and physical activity.In the present study, six shared SNPs between GLM and MLM approaches were identified on chromosomes 1 (147.6 and 167.7 Mb), 3 (29.5 and 40.4 Mb), 13 (13.2Mb), and 27 (1.9Mb) associated with body weight at hatching, a trait that is strongly influenced by maternal effects 71 .Reference 77 identified a chicken QTL on chromosome 27 in the same region as the one discovered in the present study.This region on chromosome 27 contains two genes, LIM and SH3 protein 1 (LASP1) and phosphatidylinositol-5-phosphate 4-kinase, type II, beta (PIP4K2B), which regulate the decrease in mice body weight 78,79 .Associations between gene expression at the initial ages have also been reported in chickens.Reference 80 identified a genomic region (169.8-175.3Mb) on chromosome 1 associated with body weight traits.In a study of genetic dissection of growth traits in a Chinese indigenous breed with a commercial broiler cross 81 , the major QTL of body weight was mapped to the end of chromosome 1 (173.7 Mb).Here, Vol:.( 1234567890 www.nature.com/scientificreports/associations underlying early growth stages were located on chromosome 1 between 168.74 and 172.99 Mb as indicated by both GLM and MLM analyses.Retinoblastoma 1 (RB1) was determined to be the significant gene in this position using a haplotype approach 82 .Maximum growth rate and body weight at sexual maturity and later ages are similar traits that describe mature body weight and have a high genetic correlation of nearly 1 83 .The sexual maturity in quail is between 6 and 7 weeks of age, while that in chickens is between 16 and 24 weeks of age, depending on the breed 84 .We considered quail's SNPs associated with body weight around the age of the birds in our study for comparison with the information obtained from chicken studies.The genetic regulation of body weight at maturity and later ages is a multifaceted process involving a network of genes that interact with various physiological systems, hormones, and environmental factors.The complexity of these interactions contributes to the expectation that numerous genes are involved in shaping body weight traits in quails.In a GWAS of carcass traits of Jinghai Yellow chickens 47 , five SNPs were detected for foot weight on chromosome 4 located between 75.54 and 75.67 Mb; thus, indicating that this region is located within 2.08-2.38Mb away from genes, such as family with sequence similarity 184 member B (FAM184B), quinoid dihydropteridine reductase (QDPR), and LIM-domain binding factor 2 (LDB2).In addition to these genes, Ref. 40 indicated that the F box and leucine-rich repeat protein 5 (FBXL5) has a significant influence on chicken growth traits and important biological functions.They identified this gene to be located on chromosome 4 from 72.9 to 77.9 Mb and to be strongly associated with body weight for weeks 6-16.QDPR and LDB2 were associated with shank circumference 32 and body weight in Beijing You chickens 38 .Reference 85 detected two SNPs that were associated with the body weight at the time of oviposition and are located separately at 78.8 Mb and 78.9 Mb on chromosome 4 within the FAM184B gene.The non-SMC condensing 1 complex, subunit G (NCAPG) gene, which is located within 0.1 kb downstream of the FAM184B gene, is a well-known candidate gene for body frame and carcase traits in cattle, and is considered to modulate the body and carcase weight of various breeds 86,87 .At 73.9 and 76.1 Mb, two SNPs associated with ESS 2 and ESS 1 on chromosome 4, respectively, were identified by both GLM and MLM analyses.Since ESS is positively correlated with body weight 84 , the influence of these genes on body weight is of considerable importance.Moreover, Ref. 88 reported the importance of the myostatin (MSTN) gene located on chromosome 7 on chicken body weight at 112 days of age.Similarly, we detected numerous associations underlying body weight at maturity and later ages on chromosome 7 at 22.65 Mb in both GLM and MLM analyses.The detection of QTL at later ages has also been reported in chickens 89 .The growth hormone-releasing hormone receptor (GHRHR) gene was mapped at position 1.7 cM and polymorphism analysis of this gene revealed three SNPs in the promoter region, which had a significant effect on body weight at 7, 9, 11, 13, and 17 weeks of age 90 .These findings are in line with our SNP association results.In MLM analysis, we identified 10 SNPs that explain over 10% of the phenotypic variance.This underscores their genetic importance and supports our SNP association findings.However, further investigation is needed to uncover the functional roles of these genes, shedding light on their contribution to the observed phenotypic variations at 10-16 weeks of age.
In addition to body weight, egg quality and production traits are the major selection criteria for poultry breeding.The findings of the association study on the TLE and AFE in this study hold significant implications for egg-type poultry breeding.The results have identified specific associations associated with variations in egg production traits, offering valuable insights into the genetic basis of egg production.However, the results of egg production may vary across different studies, and this variation may be attributed to the diverse durations of the egg collection period.In a GWAS investigating egg production and quality traits in chickens, Ref. 91 found the most significant SNP associated with egg number located within intron 12 of the growth factor receptor-bound protein 14 (GRB14) gene.GRB14 mRNA is highly expressed in the ovary, liver, kidney, and skeletal muscles of humans and mammals 92,93 .Similarly, our MLM suggested that an SNP associated with TLE is located on the same chromosome and position (chromosome 7; 21.46 Mb).It has been speculated that GRB14 influences egg production in layers 91 .The age of first egg is an important indicator of sexual maturation in female birds and is influenced by genetic and environmental factors 94 .Here, we detected 7 associations that might influence AFE, and they were located on chromosomes 1-4 and Z based on the GLM, and 18 associations positioned on chromosomes 1-3, 5, 8, 12, 19, and 21 based on the MLM analysis.A SNP in intron 2 of the odd oz/ten-m homolog 2 (ODZ2) gene has been previously identified to be significantly associated with chicken AFE 91 .This gene is expressed in the developing chicken brain and may affect the sexual maturity 95 .Additionally 96 , demonstrated that the gonadotropin-releasing hormone I (GnRH-I) gene is associated with chicken AFE.GnRH stimulates the synthesis and secretion of gonadotrophins, which induce steroidogenesis in the gonads, culminating in ovarian follicle growth and ovulation for egg production 97 .Reference 98 identified two novel forkhead box L2 (FOXL2) and growth differentiation factor-9 (GDF9) genes associated with egg production in the Chinese Dagu chickens.These polymorphisms play a critical role in the regulation of ovarian development in hens.This result supports our findings in the present study, in which 53 associations were detected for EPR by both the GLM and MLM analyses.
Our association study focused on external egg quality traits.Eggshell quality is a major concern in poultry breeding due to its implications for both reproductive performance and human consumption.In this study, we successfully identified several associations associated with external egg traits in both egg production stages.These findings shed light on the genetic basis of egg quality and provide valuable information for targeted breeding efforts aimed at enhancing egg production and ensuring high-quality eggs for consumers.A number of QTL regions for egg weight, egg length, and eggshell characteristics have been already identified 17 .In the present study, 49 associations were identified for EW 1,2 by both the GLM and MLM analyses.The associations for EW on chromosome 4 between 35.49 and 65.59 Mb are in line with the SNP for the same trait and chromosome at 49.28 Mb located in the shroom family member 3 (SHROOM3) gene 99 .Shroom3 is a cytoskeletal protein involved in regulating cell shape (arrangement and remodelling) in certain tissues 100 .Among the external egg traits, eggshell characteristics are the most important and play a significant role in the reproductive performance and human consumption.Reference 91 101 .Polymorphisms of eggshell organic matrix genes were considered to be related to eggshell thickness, eggshell strength, and dynamic stiffness 102 .SCNN1 has four family members, SCNN1a, SCNN1b, SCNN1g, and SCNN1d, which affect eggshell formation 103 .A study investigating eggshell quality traits identified SCNN1a, SCNN1b, SCNN1d, and SCNN1g genes on chromosomes 1 (80.03-80.04Mb), 14 (7.002-7.01Mb), 21 (2.43-2.44Mb), and 14 (7.01-7.02Mb), respectively, and associated them with egg weight, eggshell weight, eggshell percentage, eggshell strength, and eggshell thickness 101 .Our MLM analysis suggested associations for ESS 2 (chromosome 1, 78.06 Mb), ESW 1 (chromosome 14, 8.36 Mb), and EW 1 (chromosome 21, 2.58 Mb), located near the same chromosome, as stated by Ref. 101 .In addition, an association study identified three candidate genes, phosphatidylinositol-4-phosphate 3-kinase catalytic subunit type 2 gamma (PIK3C2G), inositol 1,4,5-trisphosphate receptor type 2 (ITPR2), and non-SMC condensin I complex subunit G (NCAPG) to be implicated in the dynamic eggshell quality and located on chromosomes 1 and 4 from 57.3 to 71.4 Mb 104 .PIK3C2G possesses the C2 domain and acts as lipid binding motif, ITPR2 has been shown to be important to the process of eggshell calcification, and NCAPG gene was discovered to be associated with eggshell weight for young hens in a genome-wide association analysis 104 .Moreover, one association analysis determined that another lipid-related gene, low-density lipoprotein receptor-related protein 8 (LRP8), a new member of the egg shell matrix protein family, was significantly associated with eggshell traits 105 .Our MLM findings are in support of the aforementioned studies, in which 6 associations were identified on chromosome 1 between 60.20 and 67.94 Mb for EW 2 and 5 associations were detected for EW 1 , ESS 2 , and ESW 1 on chromosome 4 between 58.94 and 69.85 Mb.
In our egg quality association study, we extended our investigation to internal aspects, such as yolk weight, diameter, color, and albumen weight.Understanding the genetic associations influencing internal egg quality is crucial due to the shift towards egg products in consumption.Through genetic marker analysis, we identified significant associations in internal egg characteristics, providing valuable insights for targeted breeding to enhance yolk weight and meet consumer preferences in the egg industry.Reference 91 identified the ataxia telangiectasia mutated (ATM) gene located on chromosome 1 (184.63Mb) to be associated with the YW trait.In the present study, based on MLM analysis, we introduced 17 associations for YW 1,2 on chromosome 1 between 67.94 and 172.99 Mb.Cathepsin D (CTSD) is another candidate gene that might affect egg characteristics 106 .Yolk formation involves cholesterol uptake and transport mediated by the very-low-density lipoprotein receptor on the membrane, and CTSD is the key enzyme regulating this process 107 .In this study, the detection of a wide range of associations (354) on all chromosomes (except for 22, 24, and 28) suggests that CTSD is a candidate gene associated with egg yolk quality traits.The albumen makes up approximately two-thirds of an egg's weight.Numerous QTLs have been reported for egg white characteristics according to the chicken QTL database (https:// www.anima lgeno me.org/ cgi-bin/ QTLdb/ GG/ index).Two promising genes, dopamine receptor D1 (DRD1) and msh homeobox 2 (MSX2) were found to be associated with albumen characteristics 108 .DRD1 encodes the D1 subtype of the dopamine receptor 109 .In birds, dopamine has been manifested to be involved in both stimulating and inhibiting prolactin (PRL) secretion, which has been illustrated to play an important role in the onset and maintenance of incubation behaviour 110 .Moreover, it has been shown to be associated with egg production traits 109 .MSX2 is a member of the msh homeobox family and is present in various embryonic tissues 108 .In the developing chick, MSX2 is expressed in the apical ectodermal ridge and the ectoderm of the genital tubercle; thus, playing a crucial role in the growth and patterning of the limb mesoderm 111 .Since MSX2 and DRD1 are involved in embryonic and ovarian development 108 , they could be treated as candidate genes associated with egg quality.Here, 23 associations were discovered for AW 1,2 located on chromosomes 1-3, 8-9, 14, 21, 23, and Z in the GLM and MLM analyses.Yolk colours are important internal egg traits that are influenced by genetic and environmental factors, most often through the components which are included in the quail's diet 112 .A total of 231 associations were found in the present study for SNPs that might influence yolk colours (L*, a*, and b* values) in both the GLM and MLM analyses.These associations were located on all chromosomes except for chromosomes 13, 22, 24, and 28.Genes underlying yolk colours have been previously reported in Galliformes.Reference 113 introduced three candidate genes, glucagon-like peptide 1 receptor (GLP1R), calcitonin receptor (CALCR), and tissue factor pathway inhibitor 2 (TFPI2), and suggested they are significantly associated with yolk colour.GLP1R is a member of the glucagon receptor family of G protein-coupled receptors found on beta cells in the pancreas and on neurons, and is involved in controlling blood sugar levels 114 .CALCR encodes a calcitonin receptor whose activity is mediated by G proteins which activate adenylyl-cyclase and is involved in calcium homeostasis, bone formation and metabolism, as well as lipid metabolism 115 .TFPI2 is a Kunitz-type serine protease inhibitor that inhibits matrix metalloproteinase activation and extracellular matrix degradation 116 .These genes are located on chromosome 5 within the 21.06-21.36Mb genomic region 113 .In the present study, based on the MLM analysis, the nearest located SNP associations were in the 18.53-25.60Mb range on chromosome 5 and responsible for YC-L* 1 and YC-b* 2 (Supplementary File S3).This similarity in positional candidate genes may indicate that our associations control yolk colour characteristics.
In this study, we applied the Bonferroni-corrected GLM and MLM to investigate the genetic basis of targeted traits.Our R 2 values conform to established statistical practices, considering multiple comparisons.The variability in R 2 values can be attributed to trait complexity and genetic diversity.Comparative analysis with existing research supports the significance of our findings, contributing to the understanding of targeted traits.
By the assessment of 44 phenotypic traits using GLM and MLM approaches, we showed the existence of SNP associations with all targeted traits, with the exception of EW 1 , AW 1 , and YC-b* 2 in the GLM and ELA 2 in the MLM analysis.The observed differences in GLM and MLM results may be influenced by factors like population structure, handling of missing data, genetic architecture, and assumptions made by each model.In agreement with previous studies, we hypothesize that multiple positional candidate genes are associated with body weight and egg-related traits.Future studies focusing on genotyping the SNPs within these candidate genes or linked markers would be of great use to improve the associations at the gene level.
To the best of our knowledge, this is the first study to perform an association analysis of the production traits of Japanese quail.A total of 383 SNPs and 1151 associations were obtained following GLM, and 734 SNPs and 1442 associations were obtained following MLM analyses.The identified associations were significantly correlated with body weight and egg production.Moreover, 35 annotated genes were introduced as candidate genes for the targeted traits based on their nearest positions.Identifying associations and candidate genes could contribute to better understanding the genetic factors controlling body weight, egg production, and quality traits in quails particularly, and chickens generally, and may accelerate the genetic progress in breeding strategies.

Figure 1 .
Figure 1.The Manhattan plot shows the association of SNP markers with ESS 2 in the GLM.Each dot represents an SNP.The figure illustrates the level of statistical significance (y-axis) as measured by the negative log of the corresponding p-value for each SNP.Each SNPs type is indicated by dots of different colours, which are arranged by chromosomal location (x-axis).The horizontal red line indicates the threshold of 5% Bonferroni genome-wide significance (p = − log 10 (5e−08)), and the underlined blue line presents a genome-wide suggestive at (p = − log 10 (0.3e−05)).

Figure 2 .
Figure 2. Quantile-quantile plot for ESS 2 trait based on the GLM analysis.

Figure 3 .
Figure 3.The Manhattan plot shows the association of SNP markers with body weight at WK 1 in the MLM.Each dot represents an SNP.The figure illustrates the level of statistical significance (y-axis) as measured by the negative log of the corresponding p-value for each SNP.Each SNPs type is indicated by dots of different colours, which are arranged by chromosomal location (x-axis).The horizontal red line indicates the threshold of 5% Bonferroni genome-wide significance (p = − log 10 (5e−08)), and the underlined blue line presents a genomewide suggestive at (p = − log 10 (0.3e−05)).

Table 1 .
117y weight (mean ± standard error) of LS and NS female birds and their F 1 and F 2 hybrids from 0 to 16 weeks of age.a-eMeans with different superscript letters are significantly different in each week of age (Tukey's HSD test, P < 0.05).WK 0-16, body weight at weeks of age.The data of LS, NS, F 1 (LS♂ ✕ NS♀), and F 1 (LS♀✕NS♂) are the same as those in our previous paper117.

Table 2 .
117y weight (mean ± standard error) of LS and NS male birds and their F 1 and F 2 hybrids from 0 -16 weeks of age.a-eMeans with different superscript letters are significantly different in each week of age (Tukey's HSD test, P < 0.05).WK 0-16, body weight at weeks of age.The data of LS, NS, F 1 (LS♂ ✕ NS♀), and F 1 (LS♀ ✕ NS♂) are the same as those in our previous paper117.

Table 5 .
Summary of chromosomes contained SNPs associations in production traits in MLM approach.Trait abbreviations are presented in Table4.Chr.no.

Table 6 .
SNP markers found in both GLM and MLM approaches.Trait abbreviations are presented in Table4; significant P values are indicated at level of **1% and *5%.